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Abstract. High-resolution R~50 000 long-slit spectroscopy of the inner knots of the highly symmetrical protostellar outflow 
HH212 was obtained in the 1-0 S(l) line of H2 at 2.12yum with a spatial resolution of ~0'.'45. At the resulting velocity reso- 
lution of ~ 6 km s , multiple slit oriented observations of the northern first knot NK1 clearly show double-peaked line profiles 
consistent with either a radiative bow shock or dual (forward and reverse) shocks. In contrast, the velocity distribution of the 
southern first knot SKI remains single-peaked, suggesting a significantly lower jet velocity and possibly a different density 
variation in the jet pulses in the southern flow compared to the northern flow. Comparison with a semi-empirical analytical 
model of bow shock emission allows us to constrain parameters such as the bow inclination to the line of sight, the bow shock 
and jet velocities for each flow. Although a few features are not reproduced by this model, it confirms the presence of several 
dynamical and kinematical asymmetries between opposite sides of the HH212 bipolar jet. The position-velocity diagrams of 
both knots exhibit complex dynamics that are broadly consistent with emission from a bow shock and/or jet shock, which 
does not exclude jet rotation, although a clear signature of jet rotation in HH 212 is missing. Alternative interpretations of the 
variation of radial velocity across these knots, such as a variation in the jet orientation, as well as for the velocity asymmetries 
between the flows, are also considered. The presence of a correlation between flow velocity and collimation in each flow is 
suggested. 

Key words. ISM : Herbig-Haro objects - ISM : individual (HH 212) - ISM : jets and outflows - stars : formation - techniques: 
spectroscopic 



r "~j 1 . Introduction portant observational constraints on these processes. In partic- 
rS ' ular, determination of line profiles and their spatial variation in 
g i Despite significant progress in our understanding of outflow ac- shocke d molecular and/or atomic gas through high-spectral res- 
tivity at the earliest stages of stellar evolution, several unsolved olution observations m the near-infrared is a powerful tool for 
problems persist, one of the most fundamental of which is the characterizing jet properties, given some reasonable assump- 
mechanism that launches and collimates protostellar jets (see tionfj about the shock physics (e g Yu et al ^pm . Despite 
reviews by Hartigan et al.EEBI Reipurth & Bally I2EU Bally, that fact> only a relatively small number of studies have used 
Reipurth & Davis |2007] Ray et al. |2003. The kinematics of high enough resolution t0 derive line profiles (e g . Zinn ecker et 
knots and bow shocks in protostellar outflows can provide im- al rj989] Carr[1993 Davis & Smith[l996] Schwartz & Greene 
[19991 Davis et al. 120011 Schwartz & Greene |2003]l. 
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Bally & Reipurth EOOTl Andrews et aL 12004) . A series of pos- 
sible explanations for the asymmetry in brightness of irradi- 
ated jets based on kinematical asymmetries has been suggested 
(Reipurth & Ballv l2"00Tl Bally & Reipurfh |20"02l They include 
the three possible combinations of equal/unequal jet velocities 
and transverse jet spread in the two arms. 

Several studies have shown the existence of a small radial 
velocity gradient across jets originating from a few Classical 
T-Tauri stars, including DGTau, RWAur, Th28, CWTau 
(Bacciotti et al. [20021 Coffey et al. [20041 [20071 Woitas et 
al. (2005 ), as well as that from a Class I protostar (HH26, 
Chrysostomou et al. 120081 ) and a Class protostar (HH212, 
Davis et al. 120001 1. As these transverse velocity gradients were 
measured consistently in most targets and with the same sign 
in both bipolar jet lobes, they were considered as observa- 
tional evidence of jet rotation. In addition, the implied rota- 
tional motions are in agreement with predictions from mod- 
els of magneto-centrifugally launched jets (e.g. Pesenti et al. 
2004), further supporting this interpretation. The search for 
such evidence is of particular importance for our understand- 
ing of star formation since protostellar outflows are thought to 
be a major channel for the removal of the excess of angular 
momentum brought into the protostellar accretion region by in- 
falling material and/or an accretion disk. The detection of such 
rotation signatures has been recently questioned in the case of 
RW Aur where it has been suggested that the disk is rotating in 
the opposite sense with respect to the bipolar jet (Cabrit et al. 
2006 ). This casts some doubts on the interpretation of the ob- 
servations in the above studies and calls for both more detailed 
observations of jet kinematics and alternate explanations for 
the transverse velocity shifts observed in jets. Alternate inter- 
pretations of transverse velocity gradients include asymmetric 
shocking, the interaction with a warped disk (Soker 2003) or 
jet precession (Cerqueira et al. l2006b . 

HH212 is one of the most remarkable protostellar out- 
flows known to date (Zinnecker, McCaughrean & Rayner 
1998). Located in the L 1630 molecular cloud at a distance 
of ~450pc, it is driven by the deeply embedded, Class low- 
mass protostar IRAS 05413-0104. Due to its clear, symmetric 
structure, HH 212 is considered to be a textbook case of a pro- 
tostellar jet. The structure of its two lobes shows pronounced 
symmetry about the driving source. Close to the source, the 
jet beams are marked by a series of bright knots with sizes of 
about 1-2 arcseconds (~ 400-800 AU) and spaced by a few arc- 
seconds. Further out, the jet is bracketed by pairs of succes- 
sively larger bow shocks (see e.g. McCaughrean et al. 2002). 
Not detected in the optical, HH212 is seen at infrared (H2 ro- 
vibrational and [Fell] lines) and (sub)millimetre wavelengths 
(e.g., CO rotational lines and SiO lines). A relatively extensive 
summary of previous observations of HH212 can be found in 
Smith, O'Connell & Davis (T20071) . 

A particularly interesting feature of HH 212 is that it is one 
of the few protostellar outflows showing some evidence of jet 
rotation, and notably the first case reported. Transverse radial 
velocity gradients of a few kms -1 have been detected in the 
H2 line at 2.12 fim from several shock knots using long-slit 
spectroscopy (Davis et al. 120001 ), although with relatively poor 
velocity resolution (~20kms~ 1 ) and angular resolution (079 
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Fig. 1. Overview of the inner part of the protostellar out- 
flow HH212 in the 1-0 S(l) line of H 2 (left, image from 
McCaughrean et al. 120021) . and blow-ups of the bow shocks 
NKl (up) and SKI (down) with the position of the slits over- 
laid. The origin of the coordinates is the approximated locus of 
the knot apex. 

pixel size). Similarly, a flattened disk-like NH3 core was de- 
tected coincident with the exciting source and with a velocity 
gradient whose direction agrees with those detected in some of 
the knots (Wiseman et al. 120011 ). The same flattened envelope 
was also observed rotating around the source in the same direc- 
tion in CO and HCO + (Lee et al. |2006ll2007l ). However, higher 
quality data are required in order to properly test the interpre- 
tation of jet rotation. 

The purpose of the present study is to probe the symme- 
try in jet velocity between the two arms of the HH212 outflow 
using the velocity structure of the inner knots as a proxy and 
test the possibility of jet rotation in HH212 from the variation 
of transverse velocities across these knots. The observations 
and their results are described in Sect.[2]and[3] respectively. In 
Sect.|4]we compare the observations with a simple bow shock 
model. Finally, in Sect.[5jwe discuss the results and implica- 
tions in terms of knot velocity structure, jet velocity asymme- 
tries, and jet rotation. 

2. Observations and data reduction 

High resolution R~50 000 long-slit spectroscopy of the HH 212 
jet bow shocks was obtained at the Gemini South Observatory 
in the nights 2003 January 12 and 13, using the near-infrared 
long-slit spectrograph PHOENIX (Hinkle et al. l2003l . This in- 
strument uses a 512x1024 InSb array, with a pixel scale of 
0.085 arcsec/pixel. The 4 pixel-wide (0734 with 07085 pixel 
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Table 1. Log of the observations. 
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Date 


Knots 


Slit positions" 


Total integration time 
(min) 


2003 Jan 12 


NK1 


perp3 


60 




SKI 


perp 


60 


2003 Jan 13 


NK1.NK2, NK4 


par 


16.7 




NK1 


perpl 


30 




NK1 


perp2 


26.7 




NK1 


perp4 


26.7 




SKI, SK2, SK4 


par 


15 



" Designation of the slit orientation with respect to the flow axis, "par" stands for parallel, "perp" for perpendicular, see Fig.fT] 



size) slit provides a nominal resolving power of 5.1 km s at 
2.12 /urn. The spectral resolution achieved as measured from 
unresolved terrestrial absorption lines in the calibration star 
spectra was ~6kms~ I FWHM, i.e. close to the nominal res- 
olution. During the two nights the seeing was on average 0745, 
as determined from the Gaussian fit to a star image profile in the 
field. This corresponds to ~ 200 AU at the distance of HH 212. 

The image presented in Fig. 1 gives an overview of the in- 
nermost knots of the HH212 jet in the 1-0 S(l) line of H2 at 
2.12jum, with the slit positions indicated and labeled. Details 
of each knot including the slit orientations and the total inte- 
gration time per orientation are given in TableQ] The two bright 
inner knots NK1 (northern) and SKI (southern) have been ob- 
served with slit positions both along and orthogonal to the jet 
axis. In particular, four adjacent slit positions with an orienta- 
tion perpendicular to the flow axis were obtained for NK1. In 
addition to NK1 and SKI, the slit parallel to the flow axis in- 
cludes as well the next two knots on both sides of the jet, i.e. 
NK2, NK4 and SK2, SK4. Note that SK3 is very faint and thus 
not visible in these data, and that there is no matching knot 
NK3 (Zinnecker et al. [T99"8l . 

Data reduction was performed as follow. Each frame was 
sky subtracted pairwise and flat-fielded. The resulting 2D spec- 
tra were corrected for bad-pixels, registered and co-added. The 
orientation of the 2D spectra with respect to the array were 
rectified row by row using the IRAF task imlintran. The an- 
gle of rotation applied, as measured from telluric calibration 
stars, was 2.72+0.05° and 2.68+0.02° for the January 12 and 
13 nights, respectively. A small bending of the 2D spectra of up 
to a few pixels was also corrected for by using the OH lines and 
a low-order polynomial fit to the row-by-row location of these 
lines given by Gaussian fit. The wavelength calibration of the 
2D spectra has been done using three bright OH sky lines in the 
raw data frames (at 2.1176557, 2.1232424 and 2.1249.592/im 
in vacuum wavelengths, from Rousselot et al.|2000). A 6th or- 
der polynomial fit was then performed for each slit position to 
determine the final wavelength solution, with residuals below 
0.03 A. Standard stars were observed in order to check the tel- 
luric absorption features, and they were found to be negligible 
in the spectral range considered. 



3. Results 

The Position- Velocity (PV) diagrams are shown in Fig. [2] 
and [3] All data are presented in the local standard of rest 
(LSR) velocity reference frame. The correction was computed 
with the IRAF task rvcorrect, assuming a conversion of 
the 2. 12 183 /mi vacuum wavelength to an air wavelength of 
2. 12125 jum. Both the high velocity resolution and the high spa- 
tial resolution of the data allow a detailed view of the velocity 
field of the H2 emission in the inner knots of HH212. In par- 
ticular, the PV diagrams with the slit oriented perpendicular to 
the flow axis have a very distinctive shape, especially for the 
NK1 knot at its widest extension (NKl-perp3 and 4 in Fig.[2j. 
It shows most of the emission at low velocity but with a tail 
toward high velocity. Also, NKl-perp3 and 4 have a parabolic 
shape towards positive velocities and a tail of triangular shape 
towards negative velocities while SKl-perp appears to be the 
reverse. This is a hint that the structure of these knots may be 
more complex than that of simple bow shocks. 

From the PV diagram of the southern flow with the slit ori- 
ented along the jet axis (Fig. [3] left panel), it is evident that 
there is a gradual redward increase in the mean knot velocities 
from knots SKI to SK4 in agreement with the data shown by 
Davis et al. (2000). Although a similar but blueshifted trend 
may exist in the northern flow, it is not immediately apparent 
in these data (Fig. 12 upper panel). The systematic increase in 
the mean radial velocity of the knots with projected distance 
from the source could be true velocity variation if a similar 
proper motion trend were observed and would therefore im- 
ply jet acceleration. Alternatively, this apparent acceleration 
could also be an indication for a bending of the flows. To dis- 
tinguish between these two possibilities will require measure- 
ment of the knots' proper motions in HH 212 (McCaughrean et 
al., in preparation). The emission between knots seen in deep 
H2 imaging (McCaughrean et al. 2002, FiglTJ is detected and 
shows a velocity dispersion of ~ lOkms^ 1 . Although the na- 
ture of this inter-knot emission is unclear, it is probably due to 
oblique shocks in the jet itself. Its increasing blueshifted mean 
velocity in the northern flow, especially apparent between knots 
NK2 and NK4, is reminiscent of the sawtooth pattern seen 
in numerical simulation of pulsed jets (e.g. Stone & Norman 
119931 1. A similar pattern can be seen between the southern knots 
SKI and SK2 but as redshifted emission. 
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Fig. 2. PV diagrams of knots of the HH212 northern flow. Upper panel: NK4-NK2-NK1 with the slit parallel to the outflow 
axis (denoted e.g. NKl-par). Lower panel : NK1 with the slit perpendicular to the outflow axis at four positions along the axis, 
denoted NKl-perpl, NKl-perp2, NKl-perp3 and NKl-perp4 from near the apex to upstream. The values and (spatial-spectral) 
locus of the mean velocities are indicated. The spatial coordinates for the slits perpendicular to the outflow are oriented positively 
towards East, i.e. essentially opposite to the NK1 brightness map of Fig.Q] 
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Fig. 3. PV diagrams of knots of the HH212 southern flow. Left : SK1-SK2-SK4 with the slit parallel to the outflow axis (denoted 
e.g. SKl-par). Right: SKI with the slit perpendicular to the outflow axis, denoted SKl-perp. The values and (spatial-spectral) 
locus of the mean velocities are indicated. Same orientation of the spatial coordinates for the slit perpendicular to the outflow as 
Fig. 12 i.e. essentially opposite to the SKI brightness map of Fig. Q] 



Unlike the case of SKI, the line profiles of NK1 are clearly 
double-peaked, as is expected for a spatially-resolved bow 
shock structure except at the bow apex and essentially inde- 
pendent of its orientation to the line of sight (e.g. Hartigan et al. 
[19901 Davis & Smith [l996l l. This is readily apparent in NKl- 
par and NKl-perp2. A double-peaked line profile can be un- 
derstood as the emission from opposite sides of the bow shock, 
with redshifted emission from receding gas from the far side 
of the bow shock and blueshifted emission from the side near- 



est the observer. The fact that the PV diagrams of SKI remain 
single-peaked at our spectral resolution suggests that the bow 
shock velocity for that knot is significantly lower than for NK1 . 
However, both NK1 and SKI have a velocity extension with 
full width at zero intensity (FWZI) of ~ 50 km s _1 which would 
indicate similar bow shock velocities, therefore suggesting a 
more complex picture. In addition, the double-peaked line pro- 
file of NK1 essentially disappears at perp3 and perp4. 
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Another possibility is that the double-peaked velocity pro- 
file of NK1 is due to the combination of two single-peaked 
shocks, namely the bow (or forward) and jet (or reverse) shock. 
This is expected in the leading bow shock of jets but also in the 
internal working surfaces (or knots) produced by pulsed jets 
where the initial velocity variation in the jet steepens into a 
sawtooth structure and both shocks develop (see e.g. Suttner et 
al. |19971 l. If the high-velocity component in NK1 is originating 
from the bow shock, as would be expected in the case of a rel- 
atively dense jet pulse, i.e. high density ratio between the fast 
and slow portion of the jet, (see e.g. Hartigan j_989 Hartigan & 
Ravmond [1993l De Gouveia Dal Pino & Benz |19941 ), it might 
explain why this component extends further ahead of the bow 
(i.e. closer to the apex) than the low-velocity component. In this 
case, the jet velocity and bow shock velocity would be compa- 
rable and the jet shock speed greatly reduced. Conversely if the 
high-velocity component is associated with the jet shock, as ex- 
pected for a relatively diffuse jet pulse, the jet velocity would 
be higher than the bow shock velocity. In such a scenario, the 
upstream part of the low-velocity component associated with 
the bow shock would be close to the radial velocity of the 
core/source with respect to the LSR, i.e. the systemic velocity 
V S yst which is + 1.6kms (Claussen et al. 119981 Wiseman et 
al.UOOD. This would imply that the northern jet lies very close 
to the plane of the sky. In any of these cases, the absence of sim- 
ilar distinct velocity components for the jet and bow shock in 
SKI would be difficult to explain with the assumption that both 
jet and counter-jet have similar velocities and density variation 
in the jet pulses. 

The differential radial velocity between the knots at the 
ends of the two oppositely directed jets (e.g. NK1/SK1) is 
~ 12kms _1 , with the northern jet approaching and the south- 
ern jet receding. For typical jet velocities (~ lOOkms -1 ), and 
with the assumption of equal jet velocities in the two flows, the 
inclination angle to the line of sight of the northern direction 
of the jet is estimated to be ~ 87°, in agreement with the 86° *l 
derived from water maser observations (Claussen et al.:1998). 
Our inclination estimate was made using the mean on-axis ve- 
locities of NK1-SK1, but the other knots give essentially the 
same result, with a maximum differential velocity with respect 
to this value of 2.0km s -1 for NK2-SK2, which corresponds 
to an angle of 86° for the same typical jet velocity as above. 
However, if one considers the systemic velocity V syst and keeps 
the assumption of equal jet velocity for the two jets, one ends 
up with a difference of ~ 5° between the two jets. For equal jet 
velocities of lOOkms 1 , the northern and southern jets should 
be inclined by 84.0° and 89.4°, respectively. Therefore, with 
the assumption that the velocity structure of the knots is exclu- 
sively due to shock gas in a bow shock structure, i.e. bow shock 
emission dominates the emission from a jet shock as is the case 
for a dense jet pulse, then the mean velocity of the knots co- 
incides with that of the jets. In such a scenario, it appears rel- 
atively clear that a significant amount of asymmetry regarding 
the orientation of propagation could exist between the two jets. 
The asymmetry can be somewhat reduced but not by much if 
one drops the assumption of equal jet velocity. This is a direct 
consequence of the fact that the mean radial velocity of knots 
are much more sensitive to variations of inclination when in- 
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Fig. 4. Sketch of the bow shock geometry and orientation. The 
observer lies in the x-z plane at an angle 9 to the z-axis. i„ and 
i, are the unit vectors respectively normal and tangential to the 
bow surface. The angle of incidence of the normal of the flow 
on to the bow surface is a. The bow shock geometry is defined 
by a function z(r). Vb is the velocity of the bow surface relative 
to the preshock material and V ps is the velocity of the preshock 
material with respect to the system rest frame. 

clinations are close to 90.0° than to variations of their (total) 
velocity. For instance, in the extreme case in which jet veloc- 
ities were of 150 and 50 km s~' for the northern and southern 
flow, the inclinations would be respectively 86.0° and 91.3°. It 
is therefore conceivable that part of the difference between the 
mean radial velocity of the northern and southern knots with 
respect to V syst could be due to a difference of jet velocities 
as well as inclination between the two arms. Note that we are 
here considering the velocity of the pre-shocked material, i.e. 
the speed at which the knots are moving as a whole, as the 
main contributor to this velocity asymmetry since for perfect 
bow shocks the constraint of equal FWZI indicates equal shock 
velocity. In the next section we consider the case in which the 
emission is only due to a forward bow shock and will check 
by bow shock modeling whether both asymmetries (orienta- 
tion and velocity) may co-exist when the velocity structure of 
the knots and not only their mean velocities are taken into ac- 
count. The short answer is yes, and the reader not interested in 
the detailed discussion can jump to the Sect. [5] The following 
section also discusses the effect of jet rotation as well as that of 
other possible mechanisms like jet precession and jet velocity 
shear on the model PV-diagrams. 

4. Semi-Empirical model of bow shock 

4.1. Model description 

We attempt to model the inner knots NK1 and SKI using a 
bow shock steady state model assuming a fixed geometry. We 
model the bow shock as the integration of individual planar 
shocks on each point of the bow surface with a shock velocity 
equal to the component of the bow shock velocity normal to 
the surface. The parameters defining the bow shock orientation 
and velocity are shown in Fig.|4] 
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Fig. 5. Best-fit HH 212 NKl PV-diagrams (perpl, perp2, perp3, perp4 and par) and knot brightness distribution without rotation 
"model A" (second row) and with rotation "model B" (third row) in comparison with the observed PV-diagrams and bow shock 
brightness distribution (upper row). The last two rows show the comparison between the observed (solid line) and modeled 
(dotted line) velocity and spatial profiles (model A). 



Our model considers a bow shock geometry of the form : 



£ - - (LX 
d~ s\d) ' 



(1) 



where s and d are constants. The case of s-2 corresponds to 
a paraboloid. Specifically, for a given bow inclination 6, a set 
of parameters {s,d) is computed by matching the projected ge- 
ometry of the bow as measured by the perpendicular slit posi- 
tion^), which are represented, e.g. in the case of NKl, by (rl, 



zl), (r2, z2), (r3, z3) and (r4, z4) corresponding to perpl, perp2, 
perp3 and perp4, respectively. The projected on-axis length of 
the bow Zmax is measured from the parallel slit, bearing in mind 
that the apex position is not well defined especially when one 
considers H2 dissociation and the inter-knot emission. The spa- 
tial dimensions of the bows NKl and SKI as measured from 
the FWHM of a Gaussian fit to the spatial profile of the spectra 
are reported in Table|2] 
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Table 2. Measured radial and axial dimensions of the NK1 and 
SKI knots. 



dimensions 


values [10 15 cm] 




NK1 








(rl.zl) 


1.87 


4.44 




(r2, z2) 


2.18 


5.92 




iri, z3) 


2.95 


8.15 




(M, z4) 


2.93 


10.0 








15 




SKI 

(nz) 

Zmax 


2.08 


3.70 

5.5 





The numerical calculation proceeds by first considering the 
bow shock at rest in a moving medium and dividing it in an- 
nuli of constant dz and radius r. The corresponding angle of 
incidence of the flow on to the bow surface is computed with 
the relation %=tan(j - a). Each annulus is divided into seg- 
ments of constant azimuthal angle if/. A single shock element 
has then a pre-shock velocity Vb (i n cos a + i, sin a)=\ \ i„ + Vy 
i t , where /„=(sin a cos if/, sin a sin if/, cos a), / f =(-cos a cos if/, - 
cos a sin if/, sin a), and Vb is the velocity of the preshock mate- 
rial impacting the bow. We assumed that after the shock V ± = 0, 
so that emission from shocked H2 arises from gas flowing along 
the bow surface with velocity Vy. This is a good approximation 
in the case of a normal jump or J-shock, as discussed in e.g. 
Schultz, Burton & Brand (2005). With this assumption, each 
individual planar element contributes to the spectrum at a radial 
velocity given by the projection of Vy onto the direction of the 
line of sight, i.e. V ra d=Vy it ■ ig, where /e=(sin 9, 0, cos 0). In ad- 
dition, in order to now consider the bow moving in a preshock 
medium itself traveling at a velocity V ps with respect to the 
system rest frame, we have to add a term - (Vb+V ps ) cos 9. 
Therefore, each segment contributes to the spectrum at a radial 
velocity : 



Vrad = Vb • ( sm a cos a cos <A sm # ~ cos a cos 6) 



V ps . COS 



(2) 



A detailed derivation of this formula can be found in 
Schultz et al. (2005) for the case in which the preshock ma- 
terial velocity in the system rest frame is zero, i.e. when the 
bow is traveling at the same velocity as the flow impacting on 
it. 

The flux contribution from an element is weighted by a 
factor proportional to its area, i.e. dA=(Rdi^)(dz/sina). The 
amount of kinetic energy from the shock that results in the 
considered molecular transition is governed by the empirical 
quantity p called the cooling function index (Schultz et al. 
2005 ), which means that the flux is proportional to V^. We 
also assume a lower shock velocity cutoff for H2 excitation at 
5 km s~' and a dissociation speed of 25 km s~' adequate for typ- 
ical molecular cloud densities (~ 10 4 cirT 3 ) and Alfven speeds 
(~ 1 kms~') (Smith |T994b . These values correspond to a mag- 
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Fig. 6. Best-fit HH 212 SKI PV-diagrams (perp and par) and 
knot brightness distribution without rotation "model A" (sec- 
ond row) and with rotation "model B" (third row) in compari- 
son with the observed PV-diagrams and bow shock brightness 
distribution (upper row). The two last rows show the compari- 
son between the observed (solid line) and modeled (dotted line) 
velocity and spatial profiles (model A). 



netic field of ~ 50 f/G. H2 emission and velocity field are com- 
puted on a grid of projected bow-axis distance z sin 6 and pro- 
jected radius R sin if/, then integrated into each slit position tak- 
ing into account its projected width. Finally, the resulting spec- 
tra are smoothed to the spatial and spectral resolution achieved 
by our measurements. 
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Considering that the shocked gas presents a rotational mo- 
tion into the x-y plane leads to the addition of a second term in 
Eq.|2]of the form : 

Vrot ■ sin if/ . sin 9, (3) 

where V mt is the velocity of rotation, assumed to be constant 
with R. Finally, we also have to include the radial velocity of 
the system with respect to the LSR. This is done by adding the 
systemic velocity V syst to the radial velocity in Eq.[2] 

Our model is empirical in the sense that no shock physics 
is included. Specifically, the emission is assumed to come from 
a thin layer which is not resolved. In other words, the cool- 
ing length of warm H2 is supposed to be smaller than the 
spatial resolution, i.e. ~3xl0 15 cm (-200AU). Smith et al. 
(2003 ) state that the approximations of a negligible cooling 
length and shock thickness hold in the case of high density. 
For HH212NK1, Tedds et al. d2002b and Smith et al. d2007b 
found a density of 10 6 and 9 x 10 4 cirT 3 in their C- and J-shock 
modeling, respectively, supporting our approximation. 

4.2. Fitting results 

Modeled PV diagrams have been computed with various bow 
shock velocity Vb, pre-shocked material velocity V ps , bow in- 
clination to the line of sight 9 and rotational velocity V rot with 
the aim of reproducing the observed PV diagrams of HH212 
NK1 and SKI separately. A least-squares fit scheme was em- 
ployed to find the set of parameters able to best reproduce those 
PV diagrams. We searched grids of parameters for the mini- 
mum of a x 2 r value. In a first step, large grids of parameters 
were searched, narrowing the parameter space to a final fine 
grid. From this first large exploration of the parameter space, it 
is evident that the outflow axes have to be almost perpendicular 
to the direction of the line of sight in order to match the separa- 
tion of the peaks in the velocity profiles for reasonable values 
of Vb. Also, since the shift of the velocity profile with respect 
to the zero LSR velocity depends both on 9 and V ps , it appears 
that 9 and/or the pre-shocked material velocity for NK1 must 
be significantly different from that for SKI. 

The fine grid was generated in the case of NK1 by varying 
Vb from 20 to 65kms~ 1 in increments of 1 kms _I , V ps from 
45 to lOOkms 1 in increments of 5kms _1 , 9 from 80° to 90° 
in increments of 1° and V rot from to 4.0km s _I in increments 
of l.Okms -1 , V rot being oriented positively counter-clockwise 
as viewed when looking north along the axis of propagation. In 
the case of SKI, the same ranges were used except for 9 from 
90° to 100°, V b and V ps which are from 15 to 35kms-' and 
from 25 to 55 km s _1 , respectively. 

Fig.|5]and|6]show the best fit PV diagrams and knot bright- 
ness distributions together with the observed for NK1 and SKI, 
respectively, as well as the corresponding velocity and spatial 
profiles. The velocity and spatial profiles are the integration of 
the PV diagrams along the spatial and velocity coordinates, re- 
spectively. The fitted parameters are reported in Table[3] The 
effect of rotation on the computed PV diagrams is shown with 
two fits, one with V rot = (model A) and one with V lot + 



Table 3. Best-fit parameters of the bow shock modeling of the 
knots NK1 and SKI. 



Parameter NK1 SKI 



Model 


A 


B 


A 


B 




w/o rot 


w/ rot 


w/o rot 


w/ rot 


en 


83 


82 


93 


92 


VUkms- 1 ) 


55 


55 


31 


31 


V ps (kms-') 


55 


50 


25 


40 


V ro , (kms- 1 ) 




2.0 




1.0 


P 


1.6 


1.7 


1.8 


1.8 


s 


1.7 


1.7 


1.9 


1.9 




28.7 


28.2 


146.7 


146.6 



(model B). The models reproduce the general shape of the ob- 
served PV diagrams, with the notable exception that the in- 
crease in velocity towards the wings of the bow shock appar- 
ent in the observed NKl-par and also SKI -par PV diagrams is 
not matched. Instead our models show the typical "spur" struc- 
ture seen in several models and observations of bow shock (see 
Section lBTTt . The blueshifted tail in NK1 and a similar red- 
shifted pattern in SKI is also not reproduced by the models. 
If present, rotation would show up as a characteristic asym- 
metric pattern in the PV diagrams, especially apparent for the 
perp3 and perp4 slit positions in NK1. Including rotation leads 
to a similar PV-diagram fit, in the chi-squared sense, as the one 
without. Therefore, we conclude that the current modeling does 
not provide a firm answer about the presence of rotation on the 
velocity level of 1-2 km s _1 . 

4.3. Transverse velocities 

In Fig.|7]we plot the barycentric (i.e. intensity weighted) trans- 
verse velocity as a function of radial distance from the jet axis 
for each perpendicular slit. The errors bars are based on the 
photon noise in the PV diagrams for which the typical peak 
signal to noise ratios are ~ 10. In order to reduce the uncer- 
tainties the barycentric transverse velocities were binned by a 
factor three. This yields an error in barycentric transverse ve- 
locities of typically 0.7 kms 1 at a radial distance of 075. The 
velocity distributions are rather flat when the slit is close to the 
cap of the bows and, in the case of NK1, become more peaked 
and redshifted towards the rear. This is not predicted by bow 
shock model A which gives a flat distribution in the case of 
no rotation, while model B leads to a trend in transverse ve- 
locity (Fig.[7]l. Note that varying the radial dependency of V ro t 
has only a marginal effect on the distribution of transverse ve- 
locity, slightly changing the slope of the trend along the bow 
axis. Rotation, if present, would therefore skew the symmetri- 
cal, sometimes peaked, transverse velocity distribution. 

Table|5]shows the slope of the transverse velocity measured 
at a radial distance of 075, defined as the difference of trans- 
verse velocity between the -075 and 075 positions. A similar 
trend can be identified for three of the slit positions in NK1, as 
well as in SKI, but only with marginal statistical significance. 
However, the fact that the amplitude of the transverse velocity 
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Fig. 8. Modeled HH 212 NKI PV-diagrams (perpl, perp2, perp3, perp4 and par) and knot brightness distribution for model C 
(second row), D (third row), and E (fourth row) in comparison with the observed (upper row). The plots in the last row show the 
comparison between the measured velocity profiles (solid line) and the modeled ones (dotted line, model E). See Tableland the 
text for a description of the different models. Note that the spatial coordinates of the PV diagrams with the slits perpendicular to 
the outflow and the abscissa of the brightness maps are oriented positively towards East. 



shifts are similar for all slit positions except perp2 is an indi- 
cation that such a trend may be real. The exception of perp2 
suggests the presence of complex possibly local fluctuations in 
radial velocities across the bow, perhaps due to turbulent mo- 
tions or density fluctuations. This could also reflect a case in 
which the poloidal velocity field is not symmetric to the jet 
axis, e.g. in the case of asymmetric bow shock wings. In other 
words, the measured slope in transverse velocities may not be 
due to rotation alone and probably includes a contribution from 
other mechanisms. 



4.4. Additional possible mechanisms 

A number of additional mechanisms could contribute along 
with the bow shock to reproduce the observed PV-diagrams 
and in particular the transverse velocity distributions of Fig. [7] 
These include 1) asymmetric bow shocks arising from a slowly 
precessing jet (e.g. Smith & Rosen 2005) or as the result of 
the orbital motion of the jet source around a binary companion 
(Fendt & Zinnecker [l998l Masciadri & Raga[2002]), 2) a high- 
velocity component which could be a Mach disk (or jet shock), 
and 3) density and/or velocity fluctuations due to turbulent mo- 
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Fig. 7. Barycentric transverse velocities as a function of radial 
distance from the jet axis measured for the different perpendic- 
ular slit positions of NK1 (top) and the single slit of SKI (bot- 
tom). Diamonds denote perpl, circles perp2, squares perp3, 
and triangles perp4. The curves correspond to model B of 
Sect. 14.21 (continuous, dotted, dashed, and dashed-dotted lines 
are for perpl, perp2, perp3, and perp4, respectively), while the 
horizontal dotted line corresponds to model A. The scale of the 
radial distance to the jet is oriented positively towards East. 

tions and/or interaction between the jet and the circumstellar 
environment (bow shock entrainment). 

In order to investigate the effect of some of these additional 
mechanisms, we modified the model developed in Section l4~T1 
We incorporated a misalignment of an angle 5 between the bow 
symmetry z-axis and the impact velocity vector Vt>. We further 
assumed that the latter rotates in the x-y-plane along the z-axis. 
Altogether, such a modification should approximate a rotating 
flow axis as this is the case for a precessing jet, but we caution 
that comparison with proper hydrodynamic calculations would 
be necessary to confirm that this is strictly equivalent to a pre- 
cessing jet. The second row of Fig. [8] shows the result of such 
a model (hereafter model C) applied to the case of NK1 with 
the best-fit parameters of model A and with 6=3°. The phase 
of rotation was also tuned to compare the observed wiggling 
pattern with a periodicity of the oscillations along the z-axis of 
2300 AU. This value has been chosen because it is the length 
of the wiggles of the inter-knot emission one can see between 
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Fig. 9. same as Fig. [7] but for model C (top), D (middle), and 
E (bottom). Diamonds denote perpl, circles perp2, squares 
perp3, and triangles perp4 observations. The curves correspond 
to models for different perpendicular slit orientations (contin- 
uous, dotted, dashed, and dashed-dotted lines are for perpl, 
perp2, perp3, and perp4, respectively). See Tableland the text 
for a description of the different models. The scale of the radial 
distance to the jet is oriented positively towards East. 



NKI and NK2 (see McCaughrean et al. l2002b . The transverse 
velocity distribution for this model is shown in Fig. [9] This 
model is able to reproduce most of the characteristics of the 
PV-diagrams and transverse velocity plots that are departing 
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from an ideal bow shock model, i.e. model A, including part 
of the east-west asymmetry of the bow shock brightness distri- 
bution. Specifically, these characteristics are the relative inten- 
sity of red- and blue-shifted peaks changing from the apex to 
the wings of the bow, i.e. from perpl to perp4, and the west- 
ward shift of the red-shifted peak in perp3 and 4 with respect 
to the bow symmetry axis. However, the blue-shifted tail seen 
in perp3 and 4 is not reproduced by such a model. In our next 
model (model D), we found that part of this feature could be 
obtained by adding a velocity shear of the jet along the x-axis 
i.e. strongly peaking on the near and far side of the bow shock 
of the form : 



Table 4. Short description of the models. 



model description 



V b (r,z) = V b 



1+0 .3.|l^ 1 + r/r4 - COS,/ ' 



(4) 



where, additionally, the velocity shear increases linearly from 
the top to the rear of the bow. The so-modeled velocity shear 
corresponds to an increase in velocity towards the border of 
the jet in the preshock region of ~ 20 % at the position of perpl 
and of ~ 50 % at the position of perp4. The results of this model 
are shown in Fig.[8}third row. Velocity shears within protostel- 
lar jets are known to exist (e.g. Solf 1987 in HH24, Bally et 
al. 120021 in HH 1/2) and are often considered as observational 
evidence of entrainment of ambient gas by the jet or its high- 
velocity component (Raymond et al. 1994] Solf 119971 Raga 
et al. 120031 1 . In addition, numerical simulations of pulsed jets 
have shown that some amount of intrinsic jet velocity shear is 
necessary in order to produce knots within the jet (or internal 
working surfaces) that have a bow shock geometry (Volker et 
al.[l999). While the jet velocity shear we introduced would be 
able to reproduce the peaking of the transverse velocity distri- 
bution seen in perp3 and 4 (see Fig.|9j, the physical mechanism 
responsible for such an increase in the velocity of preshock ma- 
terial away from the jet axis would have to be clarified. It would 
be more likely that this peaking is due to the contribution of a 
high-velocity jet shock emission. 

In Table[5] are reported the slopes of the transverse veloc- 
ity distributions of the models C and D. It is clear that neither 
of the two is able to reproduce the data slope. We then imple- 
mented a last model (model E) including a rotational velocity 
of 1 kms -1 , a value that seems to be an upper limit in order to 
approximate the data slopes. We therefore conclude that, with- 
out excluding other mechanisms able to reproduce the obser- 
vations, the model of a bow shock produced by a slowly pre- 
cessing jet with a precession angle of ~3° (or at least with vari- 
ations in the jet direction) which rotates at ~1 kms -1 (counter- 
clockwise as viewed when looking north along the axis of prop- 
agation) and additionally presents velocity shear due to jet en- 
trainment or more likely a certain amount of high-velocity jet 
shock emission would be consistent with the present NK1 data. 
The situation is less clear in the case of SKI, and multi-slit ob- 
servations similar to those used for NK1 would help to under- 
stand if the same model could also apply in this knot. 



A bow shock 

B rotating bow shock 

C bow shock + jet precession 

D bow shock + jet precession + velocity shear 

E rotating bow shock + jet precession + velocity shear 



Table 5. Transverse velocity slopes (in kms _I ) at 075 radial 
distance for the different models. Uncertainties are quoted in 
parenthesis for the data. See Tableland the text for a descrip- 
tion of the models. 



slit 


data 




model 




B 


C 


D 


E 


NK1 - perpl 


-1.3 (1.7) 


-3.3 





-0.1 


-1.6 


NK1 - perp2 


0.9(1.6) 


-3.3 








-1.6 


NK1 - perp3 


-2.0(1.4) 


-3.3 


0.2 


0.2 


-1.4 


NK1 - perp4 


-1.6(1.4) 


-3.2 


0.3 


0.4 


-1.2 


SKI - perp 


-1.0 (0.9) 


-1.6 









5. Discussion 

5.1. Knot velocity structures 

Although the observed spatial and spectral structures of the in- 
ner knots NK 1 and S K 1 of HH 212 are both qualitatively repro- 
duced by an analytical model of bow shocks assumed to arise 
from internal working surfaces in a time-dependent velocity 
(e.g. pulsating) jet, a few observed features are not predicted 
by this model. These features are (i) the increase in radial ve- 
locity dispersion from the apex towards the rear of both NK1 
and SKI knots for the slit aligned with the jet-axis and (ii) the 
lack of high- velocity tails in NK1 perp3 and perp4, as well as in 
the SKI perp PV diagrams. In the following, we will first dis- 
cuss the failures as well as possible improvements of the bow 
shock model we used and subsequently develop qualitatively 
the alternative approach mentioned in Sect.[3]which consists of 
the assumption of a dual (forward and reverse) shock. 

The simple analytical model we employed in our analy- 
sis assumed the jump conditions at the shock front (see e.g. 
Appendix in Volker et al. 1999) as has been used success- 
fully in several studies to interpret the kinematics of shocked 
H2 gas in knots (see e.g. Davis et al.[2001 ). This model is also 
consistent with the ballistic jet-driven bow shock model devel- 
oped by Ostriker et al. (2001 ) and validated with hydrodynamic 
simulations (Lee et al. 2001) which assumed conservation of 
mass flux and momentum, with shocks under strongly cool- 
ing conditions. All these models predict a "spur" structure in 
the position-velocity diagrams along the outflow axis (see also 
e.g. Hartigan et al. 1 19901 1 which contrasts with the observed 
PV-diagrams of NK1 and SKI, i.e. the increase in radial ve- 
locity dispersion from the apex towards the rear of the knots. 
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In addition, this predicted velocity structure is also in contrast 
with that observed in shocked H2 for the outer bow shocks in 
HH 212 where a similar increase in velocity dispersion towards 
the rear of the bows are seen (Davis et al. 2000). Interestingly, 
"spur" structures were observed in CO for the outermost bow 
shocks of HH 212 and both their velocity structures and spatial 
morphologies could be reproduced by such a model (Lee et al. 
|2000]Leeetai.|200l]>. 

Despite the success of pulsed jet models in reproducing 
most of the observed knot properties in HH jets, the exact na- 
ture and origin of these knots are still subject to investigation. 
While it is often observed that knots in HH jets have a clearly 
defined bow shock morphology like those in HH 111, in some 
cases (e.g. in HH 30) the knots observed along the jet axis are 
merely elongated blobs. Recently, it has been suggested that 
such knot structure could be reproduced when a large mag- 
netic field strength is included in the pulsating jet model (de 
Colle & Raga 2006). An alternative scenario for the formation 
of knots in general has also been recently proposed that in- 
cludes pressure gradients between a propagating jet and a time- 
variable cocoon (Rubini et al. 2007 ). It is also worth noting that 
other mechanisms could play a role in the formation of inter- 
nal knots such as Kelvin-Helmholtz instabilities or flares due to 
magnetic reconnections between the stellar and disk magnetic 
fields (Fendt et al. 2009]). In HH212, the presence of a dou- 
ble bow shock in NK1 as suggested by Smith, O'Connell & 
Davis d2007l l is yet another possibility which may explain the 
rising radial velocity from the apex to the rear of the bow seen 
in that knot. Although the morphology of the knots NK1 and 
SKI are reminiscent of bow shock structures when observed 
at high enough angular resolution, it is not excluded that other 
processes like those mentioned above could also contribute to 
the dynamics of shocked H2 gas. 

It is also not excluded that the departures between model 
and observations are due to some lack of shock physics in the 
model. It becomes increasingly clear that the physics of shocks 
is complex and depends on a variety of parameters includ- 
ing magnetic field strength and orientation, chemistry, grain 
physics and the effect of dust on the chemistry and dynamics 
(Brand l2007l) . Here we have modeled J-shock type bow shocks, 
i.e. no magnetic field cushioning, which are consistent with 
both numerical simulations of J-shocks (e.g. Davis & Smith 
1996, Fig. 9, slit 4, 75deg.) and analytical models like those 
mentioned previously. However, line ratio analysis (Zinnecker 
et al. |1998l Tedds et al. 2002) suggests that the shocks may 
be of the C-type. Numerical simulations of C-shocks in bow 
shocks usually yield PV-diagrams which have some similarities 
with what we observed with the slit on-axis. For example, in 
Davis & Smith ( 1 19961 1 the shape of the modeled C-shock PV di- 
agrams (Fig. 10, slit 4) exhibits a less pronounced "spur" struc- 
ture like that apparent in J-shock PV diagrams. Some features 
such as the increase in the width of the gap between the two 
velocity peaks towards the wings of the bow seems to match 
our observations better. However, other features predicted by 
C-shock models like the increase in the velocity dispersion of 
each velocity peak towards the bow shock apex are not ob- 
served. Recently, Smith, O'Connell & Davis ( 120071 ) have mod- 
eled integral field (low resolution) spectroscopy observations 



between 1.5 yum and 2.5 /vm of NK1 and SKI as both C- and 
J-type bow shocks. Both the excitation and multi-wavelength 
morphology of the knots are found to be in slightly better agree- 
ment with J-type bow shocks with physical parameters (bow 
velocity and geometry) globally consistent with those we de- 
rived in this study. In addition, the existence of strong [Fell] 
emission in the inner knots of HH212 (Zinnecker et al. 1998, 
Caratti o Garatti et al. 2006 ) would be consistent with the pres- 
ence of dissociative J-type shocks at least at the cap of the 
bows. Therefore, although J-type shocks seem to be more ap- 
propriate to explain the integral field observations, the inter- 
mediate case of a bow shock with dissociative J-shocks at the 
cap and C-shocks in the wings or alternatively a J-shock with 
a magnetic precursor (Giannini et al. 2004) could also be con- 
sidered. 

As already mentioned in Sect. [3] a more likely alternative 
would be that the observed PV-diagrams in NK1 arise from 
a dual (forward and reverse) shock structure, i.e. a combina- 
tion of a jet shock (or reverse shock or Mach disk) and a bow 
shock. In such a case, there are two possibilities. First, if the 
high-velocity component is associated with the forward bow 
shock, then the velocity of the jet shock would be significantly 
lower than the bow shock speed which, under the assumption 
of ram pressure equilibrium in the working surface, would indi- 
cate that the jet density is larger than the density of the preshock 
material (e.g. Blondin et al. 1990 De GouveiaDal Pino & Benz 
1994 Hollenbach 1997). For a pulsed jet this means that the 
density of the fast jet material is higher than that of the slow 
jet material, i.e. there is no mass flux conservation in the jet 
pulses (Hartigan & Raymond 119931 De Gouveia Dal Pino & 
Benz |199"4j ). Also, in that case, the jet velocity is of the order 
of the forward bow shock speed. This scenario would be con- 
sistent with the spatial separation of the two components, with 
the higher velocity component located further ahead of the knot 
than the low velocity component. However, the decrease of ra- 
dial velocity of the high velocity component towards the apex 
would be difficult to explain in a standard bow shock analysis 
as discussed above. In addition, the forward shock would be 
expected to be stronger than the jet shock and would therefore 
likely dominate the knot emission leading to a brighter high- 
velocity component which is not observed. The second possi- 
bility would be that the high-velocity component is associated 
with the jet shock. In that case, the decrease of radial velocity 
towards the knot's apex would presumably be due to the jet de- 
celeration, and molecular dissociation near the bow apex could 
account for the predominance of low-velocity bow shock emis- 
sion in the wings. This scenario would thus imply that the jet- 
to-preshock material density ratio is significantly smaller than 
unity (under the same assumption of ram pressure balance as 
before), which implies mass flux conservation in the jet pulses. 
This would be the case if the jet density varies inversely with 
the jet velocity in the pulses (Hartigan & Raymond l 19931 Stone 
& Norman [1993]). In summary, if the knot is associated with a 
"dense pulse", we would only observe the forward shock with a 
velocity similar to that of the jet. Conversely, if the increase in 
jet velocity giving rise to the knot corresponds to a decrease of 
jet density ("diffuse pulse"), then both jet and forward shocks 
are present with similar brightness and the velocity of the jet is 
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higher than that of the jet shock which in turn is higher than the 
bow shock speed. This is the case we would favor for NK1 . In 
this context, the analysis performed in Sect.|4]would be equiv- 
alent to the case of NK1 and SKI both arising from a "dense 
pulse". 

Disentangling the scenario of a single, velocity-resolved, 
bow shock and that involving a combination of velocity unre- 
solved jet and bow shocks could in principle be achieved by 
measuring the excitation temperatures in each low and high ve- 
locity component of NK1 . As the jet shock would be expected 
to have a significantly higher temperature than the bow shock, 
differences in excitation between the two velocity components 
would favor a scenario involving a dual shock rather than a bow 
shock structure alone. As already pointed out in Sect.[3] the 
single-peaked velocity profile in SKI suggests the presence of 
velocity asymmetries and possibly also different density varia- 
tion in the jet pulses between the two arms of HH212 (see fur- 
ther discussion in Sect. 15. 2b . A further alternative would be that 
the jet shock in SKI is completely dissociative or that the jet 
is mostly atomic in nature due to previous dissociative shocks 
(see discussion in Yu et al.[2000). Further observations and de- 
tailed modeling are needed to test all these possibilities. 

5.2. Asymmetries in the jet : velocity & orientation 

Our bow shock modeling presented in Sect.|4]identified signifi- 
cant kinematical differences between the northern and southern 
inner flows : a higher velocity of the northern flow compared to 
southern flow by a factor of ~ 2, as seen in the velocity of the 
bow shocks with respect to the preshock material as well as in 
the velocity of the preshock material. Our model yields bow 
shock velocities with respect to the LSR of ~ HOkms 1 for 
NK1 and ~ 55 km s -1 for SKI corresponding to bow shock ve- 
locities with respect to the preshock material of 55 km s _I and 
31 kms -1 , respectively (Table[3|). There is therefore a factor of 
about two difference in velocity between the two lobes. In ad- 
dition, a factor of two may likely be a lower limit as the bow 
shock velocity in SKI would have to be smaller than the value 
reported above in order to reproduce the single peaked line pro- 
file. 

On the other hand, as already mentioned in Sect l5.ll if the 
observed velocity structure in NK1 is due to dual forward and 
reverse shocks, jet velocity asymmetries between the two arms 
of HH212 would also explain the apparent lack of a reverse 
shock in SKI in the case that this knot is associated with a 
"dense pulse" and NK1 with a "diffuse pulse". This scenario is 
particularly attractive as it would imply similar knot velocities 
for the northern and southern arm and would thus lead to knots 
at equidistant locations from the driving source as observed. 

A similar velocity asymmetry is also found in the kinemat- 
ics of H2 shocked gas of the outer bow shocks NB 1 and SB 1 
for which NB 1 shows a factor of ~ 2 larger velocity disper- 
sion than SB 1 and only NB 1 is resolved into a double peaked 
velocity distribution (Davis et al. 2000). If this velocity asym- 
metry is also present in the proper motion of the knots and bow 
shocks, then this has to be attributed to a true velocity asym- 
metry, otherwise the two lobes may exhibit significantly dif- 



ferent inclinations. As already discussed in Sect.[3j both asym- 
metries (velocity and orientation) could be present at the same 
time, a situation also suggested by SiO observations (Lee et 
al. [2007 1. Support for that interpretation is also provided by 
the fact that HH212 already exhibits a few degrees misalign- 
ment between opposite flows in the plane of the sky (Smith et 
al. l20071 >. However, a misalignment between the flows of a few 
degrees in the direction of the line of sight would be insufficient 
to accommodate both the asymmetry in velocity of the knots 
and bow shocks and their apparent symmetry in projected dis- 
tance from the source. This puzzling situation could perhaps be 
resolved if the flow speed would have relatively sudden varia- 
tions, but such a hypothesis seems unlikely. 

Hirth et al. d 1994b suggested that different flow speeds on 
the two opposite sides could be connected with differences 
in flow width/opening angles. High spatial resolution imag- 
ing is required to measure reliable values of opening angle 
and/or width of jets and therefore properly test such a possi- 
ble connection. From the recent literature, there are at least two 
cases for which such a trend seems to be confirmed : HH 30 
(Hartigan & Morse 12007b and RWAur (Lopez-Martin et al. 
2003, Woitas et al. 2002). In particular, there appears to be a 
correlation between flow speed and collimation : the larger the 
jet width/opening angle the faster the jet. This is also the case 
in HH212 for which the northern flow is faster and less col- 
limated than the southern one. A high-velocity 12 CO (J=2-l) 
map clearly shows that the blueshifted outflow is less colli- 
mated than the redshifted outflow (Lee et al. 2006, Fig. 6b). 
This is also seen in the 12 CO (J=3-2) map (Lee et al. [2007, 
Fig. 4). If confirmed in other outflows, such an empirical law 
would need to be explained by models. Noticeably, such a cor- 
relation between flow speed and collimation would be opposite 
to what would be expected from a jet with unequal jet velocities 
and identical transverse jet spread in the two arms (Reipurth & 
Ballv l200T1 Bally & Reipurth I2002T Different mass load could 
possibly explain such a trend, a possibility that needs to be fur- 
ther investigated. 

One would like to understand if such velocity asymme- 
tries are intrinsic, i.e. originating from the central engine, or 
propagation related. If they were due to asymmetries in den- 
sity in the ambient medium in the two lobes, then one would 
expect to observe an asymmetry in brightness, which is not ev- 
ident in HH212. There is also no evidence of differences in 
density between the two arms of the outflow from CO obser- 
vations (Lee et al. 2007). The existence of different intrinsic 
jet velocities seems thus unavoidable and the reason for this 
would need to be better understood. For example, it is con- 
ceivable that the regions of jet launching on opposite sides 
of the outflow source are located at significantly different ra- 
dial distances from the source. However, such differences in 
jet launching radius would have to be unrealistically large in 
order to accommodate a factor two in outflow speed ratio, un- 
less additional asymmetries in the mechanism responsible for 
jet launching/collimation are present (like e.g. in the ambient 
medium close to the disk which could affect the lever arm or in 
the magnetic field). In this respect, multipolar magnetic fields, 
as opposed to dipolar magnetic fields usually assumed in mod- 
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els, may be promising in order to produce such differences be- 
tween opposite sides of a jet. 

5.3. Transverse velocity gradients and jet rotation 

Our observations show the presence of transverse velocity gra- 
dients - although with limited statistical significance - with the 
same sign in the knots NK1 and SKI, and with values of 1 - 
2km s _1 at a radial distance of 075 from the jet axis (Table[5]l. 
While these transverse velocity gradients could be interpreted 
as evidence of jet rotation, with values that would agree with 
predictions from the current MHD jet models (Ferreira et al. 
2006 ), we will also discuss the possibility that other mecha- 
nisms could reproduce such gradients, including jet precession 
and a jet from a companion in orbital motion. 

A tentative jet rotation signature in HH 212 was reported by 
Davis et al. (2000 ) from the measurement of a transverse veloc- 
ity gradient in NK1 that would be consistent with a clockwise 
(when looking along the northern direction of the jet propa- 
gation) rotation speed of ~ 1.5kms _I . However the latter in- 
terpretation is not supported by several other observations. The 
radial velocity gradient measured by Davis et al. (2000]( in NK1 
is not consistent both with that measured in SKI by the same 
authors and with the velocity gradient measured in the molec- 
ular envelope around the exciting source IRAS 05413-0104 
(Wiseman et al. 120011 Lee et al. 120071 ) . In contrast, the trans- 
verse velocity gradients measured in this work have both the 
same sign in NK1 and SKI (identical to Davis et al. in SKI, op- 
posite in NK1), i.e. consistent with a counter-clockwise sense 
of rotation as viewed when looking north along the axis of 
propagation. The sign of of such transverse velocity gradients is 
consistent with the direction of the velocity gradient measured 
in the embedded source, apparently supporting the jet rotation 
scenario. 

However, while SO (sulfur monoxide) observations sug- 
gested the detection of a transverse velocity gradient of a few 
km s -1 in knot SK4 and knot SN (SO knot closer to the source 
than the H2 knot NK1) with a sign identical to what we mea- 
sured in both NK1 and SKI (Lee et al. 120061 120071 recent 
SiO observations have failed to detect velocity gradients on 
that level for the innermost SiO knots (Codella et al. 120071 ). 
In addition, our detailed PV diagrams of SKI and NK1 with 
perpendicular slits (SKl-perp, see Fig. [3] NKl-perpl, 2, 3 and 
4, see Fig.O do not show a clear signature of jet rotation as 
would be expected from our model. In the case of a bow shock 
only (or a "dense pulse"), we have shown in Sect. [4] that both 
the structure of the PV diagrams and the variation of barycen- 
tric velocities as a function of distance to the jet axis, while 
still consistent, do not provide a completely satisfactory match 
with the model that includes rotation. A similar conclusion can 
be drawn if there is a dual (forward and reverse) shock. In the 
case of the jet shock being related to the low-velocity compo- 
nent ("dense pulse"), a transverse velocity gradient would be 
present in NK1, which would have the same sign as the trans- 
verse velocity gradient seen at low velocities in SKI. However, 
the "dense pulse" scenario is not the favored case for NK1 (see 
Sect. 15. It . In addition, there would be no significant transverse 



velocity gradient in the case in which the jet shock is associated 
with the high-velocity component ("diffuse pulse") because of 
the symmetry of the kinematic al structure of this component 
with respect to the flow axis (see NKl-perp3 and 4 in Fig.[2j. 
Therefore we should be cautious in interpreting our data as the 
result of jet rotation in HH212. 

An alternative interpretation of transverse radial velocity 
shifts could be jet precession (Cerqueira et al. 2006). Thus, in 
principle jet rotation could be mimicked by jet precession, the 
reason being that the same sign of the transverse velocity gradi- 
ent would be produced for knots on opposite sides of the flow. 
However, recent hydrodynamical computations of precessing 
jets disagree with the Cerqueira et al. result (Smith & Rosen 
2007). Indeed, our own simple model described in Sect. 14.41 
disagrees, too, and would further suggests that both jet pre- 
cession and jet rotation could be operating at the same time 
in HH212. Wiggles are often observed in protostellar jets and, 
as discussed in Reipurth et al. (2000), those may be produced 
by several mechanisms. Anglada et al. (2007) have shown that 
the wiggles of the jet HH 30 would be consistent with both a 
precessing jet and a jet originating from a companion in orbital 
motion (Fendt & Zinnecker |19981 Masciadri & Raga 2002|. As 
the latter scenario would likely produce velocity gradients with 
identical directions in both jet and counterjet, it would also be 
another possible interpretation for both the transverse velocity 
shifts and the wiggles observed in HH212. However, as for jet 
precession, proper hydrodynamical computations are needed to 
confirm that this mechanism can efficiently produce velocity 
shifts in the jet knots. 

In principle, the morphology of the jet would be a way 
to distinguish between these two scenarios. A precessing jet 
would produce a point symmetric structure while a jet from 
a companion in orbital motion would be revealed by a mirror 
symmetry structure. While no such obvious structure can be 
observed in the H2 images at 2. 12/rni, a point symmetric struc- 
ture is observed in the associated CO outflow (Lee et al. 2007), 
supporting the possibility of a precessing jet in HH 212. On the 
other hand, the transverse gradient in brightness of NK1 and 
SKI (both brighter at west), corresponding to a transverse gra- 
dient in excitation as recently shown by Smith et al. ([2007), 
would be more consistent with the orbital motion scenario. 
However, such a transverse gradient would also be naturally 
explained by the transverse ram pressure produced by either a 
westward motion of the outflow source or by an eastward drift 
of the ambient cloud material. 

Clearly more detailed observations and modeling are 
needed in order to identify which one of these mechanisms is 
prevalent in HH212, and whether some of them (e.g. rotation 
and precession) are operating together. 

6. Conclusions 

In this paper, we have investigated the velocity structure of the 
inner H2 knots in the HH 212 bipolar jet. The main conclusions 
of our study are summarized as follows : 

(1) The kinematics of the knots NK1 and SKI can be in 
part reproduced quantitatively by a model of bow shock. An 
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alternative qualitative approach would be that of a dual forward 
and jet shock. The bow shock model is a particular case of this 
alternative approach when the mass flux in the jet pulse is not 
conserved ("dense pulse"). 

(2) Both the mean knot velocities and the bow shock modeling 
suggest significant asymmetries in velocity and orientation 
between the two lobes. In addition, the knot NK1 could 
be produced by a relatively "diffuse pulse" and SKI by a 
"dense pulse". We note the existence of a possible correlation 
between flow speed and collimation with the larger the jet 
width/opening angle the faster the jet. 

(3) Both the PV diagrams and the variations of radial velocity 
with distance to the jet axis are not well reproduced by a model 
of rotating bow shock alone. Although possible alternatives 
include jet precession alone and oblique shocks due to e.g. an 
s asymmetric poloidal velocity, we found that a combination 
of jet rotation and jet precession in addition to velocity shear is 
somewhat more consistent with the data. 
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